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Abstract. A state-of-the-art 3D source imaging technique is used to extract the 3D 
two-pion source function in central and mid-central Au-|-Au collisions at y/sNN = 
200 GeV. The source function indicates a previously unresolved non-Gaussian tail in 
the directions of the pion pair transverse momentum (out) and along the beam (long). 
Model comparisons give robust estimates for several characteristics of the emission 
source, including its transverse size, its mean proper breakup time r and its emission 
duration At. These estimates are incompatible with the predictions for a first order 
phase transition. However, they point to significant relative emission times which could 
result from a crossover phase transition. 



1. Introduction 

Lattice calculations indicate a rapid transition from a confined hadronic phase to a 
chirally symmetric de-confined quark gluon plasma (QGP) at the critical temperature 
Tc ~ 170 MeV [1]. Such a plasma is produced in energetic Au-|-Au collisions at the 
Relativistic Heavy Ion Collider RHIC. However, it is currently unclear whether the 
phase transition to the QGP phase is first order (AT = 0) or a rapid crossover refiecting 
an increase in the entropy density associated with the change from hadronic (dn) to 
quark and gluon (dq) degrees of freedom. 

An emitting system which undergoes a first order phase transition is predicted 
to have a large space-time extent [21 [3]. This is because the transition "softens" the 
equation of state (ie. the sound speed Cg ~ 0) in the transition region, and this delays the 
expansion and considerably prolongs the lifetime of the system. A smaller space-time 
extent has been predicted for systems which undergo a crossover transition [3]. 

To search for a prolonged lifetime, it has been a common practice to measure 
the widths (R) of the emission source function (assumed to be Gaussian) in the 
out- side- and long-direction {Rout, Rside and Riong) of the Bertsch-Pratt coordinate 
system Here, the prediction is that Rout I Rside » 1, for systems which 

undergo a first order phase transition, [2l [3]. This is illustrated in Fig. (H 
where the values obtained from a hydrodynamical model calculation [3] are plotted 
as a function of energy density (in units of TcSc-, s is the entropy density). 
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These rather large ratios (cf. Figs. [T^ 
and c) have served as a major motivating 
factor for experimental searches at several 
accelerator facilities [1]. No evidence for 
a prolonged lifetime were found by these 
studies and the reported Gaussian source 
functions are spheroidal with _Rout ~ -Rsidc 
in the longitudinally co-moving system. 

A crossover transition can be mod- 
eled by varying the width AT, of the tran- 
sition region. Figs. [1^ and b show that 
the magnitude of Rout/Rside is consider- 
ably reduced (by as much as a factor of 
four) when calculations are performed for 
AT = O.lTc. Thus, the space-time extent 
of the emitting system is significantly in- 
fluenced by the cross over transition and 
this could lead to less prominent signals 
which require more sensitive methods of 
detection. Indeed, a recent study with a 
ID source imaging technique has observed a long non-Gaussian tail in the radial source 
function and attributed it to possible lifetime effects [Ej. 

In this contribution, we report on recent efforts to study the 3D two-pion source 
function via a new state-of-the-art technique proposed by Danielewicz and Pratt [6]. 
Namely, the 3D correlation function is first decomposed into a basis of Cartesian surface- 
spherical harmonics to extract the coefficients, also called moments, of the expansion. 
They are then imaged or fitted with a trial function to extract the 3D source function, 
which is then used to probe the emission dynamics of the pion source [71 [8]. 

2. Analysis Method 

The 3D correlation function C(q) = Afgd(q)/Abkg(ci) "was obtained by taking the ratio 
of the 3D relative momentum distribution for ti^tt^ and tt^tt^ pairs in the same event 
A^fgd(q) and those from mixed events Abkg(q); where q = ^p^^p^) where pi and p2 are 
the momentum 4- vectors in the pair center of mass system (PCMS). 

The 3D correlation function C(q) was expanded in a Cartesian harmonic basis [6] 
to obtain the moments 

C(q) - 1 = i?(q) = E (1) 

l,ai...ai 

where / = 0, 1, 2, . . ., = x, ?/ or z, ^aj...^, (f^q) are Cartesian harmonic basis elements; 
(fiq is the solid angle in q space); R^^^ aS^) Cartesian correlation moments given 
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Figure 1: Rout/Rside ^ ^ function of the initial 
energy density for an expanding fireball. (a,b) 
are for dq/dn = 37/3, (c,d) for dq/dn = 3. 
The thick lines in (a,c) are for AT = 0, in (b,d) 
for AT = O.lTc. Thin lines show results for 
an ideal gas case. Solid lines show results for 
T = 0.7Tc, dotted for T = 0.9Tc, dashed for 
T = Tc. The figure is taken from Ref. [3]. 
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by Eq. 



); and q is the modulus of q. 
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Here, the coordinate axes are oriented 
so that z (long) is parallel to the 
beam direction and x (out) points in 
the direction of the total transverse 
momentum of the pair. For this analysis, 
Eq. ([1]) was truncated at / = 6 and 
expressed in terms of its 10 independent 
even moments: RI2, R^2i ^xAi 

^x2y2^ -^x6' -^y6' ^x4y2 ^x2yA i-^: 

is shorthand for Rl^ [6]); odd moments 
were checked and found to be consistent 
with zero [within statistical uncertainty] 
as required by symmetry considerations; 
higher order moments (for / > 6) were also 
checked and found to be negligible. 

The independent moments are shown 
as a function of q in Fig. [2l They were 
obtained by fitting the truncated series 
to the measured 3D correlation function 
with the moments as the parameters 
of the fit. It is noteworthy that the 
very good agreement between R^{q) and 
R{q) (shown in panel (a)) points to 
the absence of any significant angular 
acceptance issues. It also attests to 
the reliability of the moment extraction Figure 2: Experimental correlation moments 
technique. i?°(g) and i?(g) both represent R^{q) for / =0, 2, 4, 6. Panel (a) also 
angle-averaged correlation functions, but shows a comparison between R°{q) and R{q). 
R%q) is obtained from the 3D correlation Systematic errors are less than the statistical 
function via Eq. ([2]) while R{q) is 
evaluated directly from the ID correlation 
function as in Ref. [5]. 

The 3D source function S{r) is obtained from the moments via imaging or fitting. 
This is made transparent by the observation that, in analogy to Eq. ([T]), ^(r) can also 
be expanded in a Cartesian Surface-spherical harmonic basis 
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errors. The sohd curves indicate the Hump 
function Eq. ([6]) fit. 
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If the series for R{q) and S'(r) are now substituted into the 3D Koonin-Pratt equation 



C(q) - 1 = R{q) = J drK{(i,r)S{r) 



(4) 
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then, the following set of ID expressions (Eq. ^) [B] which relate the correlation 
moments Raj^,,,aiil) the source moments are obtained; 

RL..aM) = 4vr/ drr'Ki{q,r)Sl^,,,^^{r). (5) 

Here, it should be noted that S{r) gives the probability of emitting a pair of particles 
with a separation vector r in the PCMS and the 3D Kernel, K{c{, r), incorporates both 
the Coulomb force and the Bose-Einstein symmetrization. 

The ID imaging code of Brown and Danielewicz [9] was used to numerically 
invert each correlation moment i?^^ to extract the corresponding source moment 

S^^ ^^{r). The latter were then combined as in Eq. ([3]), to give the source function. 
The 3D source function was also extracted via direct fits to the 3D correlation function 
with an empirical Hump function given by 

SH(r.,r„r.) = XeM-U^ + + " ^^^^ + ^ + 1^)1' 

^'■xs ys ^^zs xl yl zl 

where A, ro, r^s, r-ys, r^s, r^i, ryi, r^i are fit parameters and /s = 1/[1 + (r/ro)^], 

/i = 1 — fs- This procedure corresponds to a simultaneous fit of the ten independent 

moments. The solid curves in Fig. [2] show the results of such a a fit. They indicate that 

the 8-parameter Hump function achieves a good fit to the data (x^/ndf=1.4). 



3. Source Image and its interpretation 



Figures [3](a)-(c) show the source function 
profiles S(r2;), S{ry) and S{r;,) obtained 
via fitting (line) and source imaging 
(squares). S(r^.) is characterized by 
a long tail, which is resolved up to 
~60 fm; S{ry) and S(r^) are resolved up 
to ~25 fm. The corresponding correlation 
profiles obtained by summation of the 
data (circle), fit (line) and image (square) 
moments are shown in Figs. [3](d)-(f). 
The broader S{rx) is associated with the 
narrower C^q^) (Fig. [3](a) and (d)), as 
expected. 

The extended tail lies along the pair 
total transverse momentum. Thus, the 
relative emission times between pions, 
as well as the source geometry, will 
contribute to S(ra;). The source lifetime 
contributes to the range of S(r^) and S(ry) 
reflects its mean transverse geometric size. 
The difference between S{r^) and S(ry) is thus driven by the combination of the emission 
time difference, freeze-out dynamics and kinematic Lorentz boost. 
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Figure 3: Source function profiles S{rx), S{ry) 
and S{rz) (left panels) and their associated 
correlation profiles C{qx), C{qy) and C(g^) 
(right panels) in the PCMS. The bands indicate 
statistical and systematic errors. 
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To aid interpretation of the source 
function, the event generator Thermina- 
tor [in] was used to predict source func- 
tions for several model scenarios. Ther- 
minator provides thermal emissions (in- 
cluding all known resonance decays) from 
a longitudinally oriented, boost invariant 
cylinder of radius pmax- A differential fluid 
element is a ring defined by cylindrical co- 
ordinates z and p; it breaks up at proper 
time r in its rest frame or at time t in the 
lab frame, where t'^ = + z'^. The freeze- 
out hypersurface is given by r = tq -|- ap, 
where Tq is the proper breakup time for 
p =0 and a represents the slope of the 
freeze-out hypersurface in p-r space; for 
a > particles at small p's are emitted 
at earlier times i.e inside-out "burning"; 
outside-in "burning" occurs for a < 0. 

Calculations were performed in Blast- 
Wave mode for a set of parameters tuned 

to fit charged pion and kaon spectra. In addition, transverse expansion [governed by 
a radial velocity Vr semi-linear in p, i.e. Vr{p) = (p/Pmax)/(p/Pmax + I't) {vt = 1.41)] 
was assumed and a was varied. The mid-rapidity pion pairs so obtained were also 
generated with the effects of all known resonance decay processes on and off. The 
resulting pairs were then transformed to the PCMS, as in the data analysis, to obtain 
S{fx,y,z) distributions for comparison with the data. 

Figure H] shows that the 3D source function generated by Therminator calculations 
(open triangles) with tq = 8.55 fm/c, pmax = 8.92 fm and a = —0.5 with resonances on, 
reproduce S{ry) but do not fully account for the long tails in x and z ie. the latter are 
longer than the Therminator source profiles. For the same parameters, the calculations 
underestimates S(ra;), S>{ry) and S(r^) when resonances are turned off (solid triangles). 
Reasonable attempts to fit the distributions by only increasing tq or with a > failed, 
suggesting substantial contribution from pion pairs with significantly longer emission 
time differences. 

An alternative approach to lengthen the distribution of time differences between 
pion pairs is to sample them from a family of hypersurfaces defined by a range of values 
of proper breakup times r'. One such parametrization consists of replacing r by r' 
chosen from an exponential distribution dN/dr' = ^^^^^^ exp[— (r' — r)/Ar], where the 
width of the distribution At represents the mean proper emission duration. Figure H] 
shows that this approach, with At = 2 fm/c (open circles), leads to a fairly good match 
to the three observed source profiles. 
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Figure 4: Source function comparison between 
Therminator calculation and image for (a) 
S(r^), (b) S{ry), (c) S{r,) in PCMS. Panel 
(d) compares AtLCM from Therminator events 
with various assumptions for At and resonance 
emission. 
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Figure 111(d) summarizes the relative emission time distribution in the LCMS, 
AtLCM, for pion pairs from events with the parameterizations indicated. For a fixed 
To = 8.55 fm/c (Ar = 0) and resonance decays excluded, the distribution AtLCM is 
narrow, (|AtLCM|) = 2.4 fm/c. The addition of resonance decays adds a long tail and 
gives (I AtLCM I) = 8.8 fm/c. Replacing r with the exponential distribution r' with 
Ar = 2 fm/c, results in a AtLCM distribution which is significantly broadened to give 
(|AtLCM|) = 11-8 fm/c. The wider distribution of time delays is needed to reproduce the 
source distributions. This implies a non-zero proper emission duration in the emission 
rest frame. 

Figure H] shows that substantial time differences AtLCM are required to account 
for the source distensions; however, the interplay between proper time and breakup 
dynamics is model dependent. Nevertheless, the picture which emerges from the 
Therminator model comparison is consistent with an expanding fireball (pmax = 8.92 fm) 
with proper breakup time Tq ~ 9 fm/c, which hadronizes and emits particles over a short 
but non-zero mean proper emission duration Ar = 2 fm/c. Such a short time duration 
is incompatible with the predictions [21 [3] for a first order phase transition. 

4. Conclusions 

In conclusion, a novel three-dimensional source imaging technique has been used to 
extract the 3D pion emission source function in the PCMS frame from Au-fAu collisions 
at a/s/v7v = 200 GeV. The source function has a much greater extent in the out (x) 
and long (z), than in the side (y) direction. Therminator model comparisons suggest 
an emission source (pmax = 8.92 fm) burning from outside in with proper lifetime 
ro ~ 9 fm/c and a mean proper emission duration Ar ~ 2 fm/c. These emission 
characteristics are incompatible with the predictions for a first order phase transition. 
However, they point to significant relative emission times ((|AtLCM|) ~ 12 fm/c, 
including those due to resonance decay) which could result from a crossover phase 
transition. 
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